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ABSTRACT 

We present FUSE observations of the line of sight to WD 0439+466 
(LS V+46 21), the central star of the old planetary nebula Sh 2-216. The FUSE 
data shows absorption by many interstellar and stellar lines, in particular D I, 
H 2 (J = - 9), HD J = - 1, and CO. Many other stellar and ISM lines are 
detected in the STIS E140M HST spectra of this sightline, which we use to de- 
termine iV(H I). We derive, for the neutral gas, D/H = (0.76 ± ^ 2 )xl0" 5 , O/H 
= (0.89 ± °;^)xl0- 4 and N/H = (3.24 ± ° : |)xl0- 5 . We argue that most of the 
gas along this sightline is associated with the planetary nebula. The low D/H 
ratio is likely the result of this gas being processed through the star (astrated) 
but not mixed with the ISM. This would be the first time that the D/H ratio has 
been measured in predominantly astrated gas. The O/H and N/H ratios derived 
here are lower than typical values measured in other planetary nebulae likely due 
to unaccounted for ionization corrections. 

Subject headings: ISM: Abundances - - ISM: Evolution - - Ultraviolet: ISM 
- Stars: Individual (WD 0439+466, LS V+46 21): ISM — Planetary Nebulae 
(Sh 2-216) 



1. INTRODUCTION 



Stars with intermediate mass (1-8 M ) evolve through the asymptotic giant branch 
(AGB) experiencing large mass-loss rates, replenishing the interstellar medium with metals, 
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on their way to become white dwarfs with M < 1.44M Q . When the stars become hot enough 
they photoionize the material ejected during the AGB phase, giving origin to a planetary 
nebula (PN). The gas in the PN reflects the stellar pro cessing that occurred du ring the stellar 
life, for species such as helium, nitrogen, and carbon (jStanghellini et al.ll2006l ). while oxygen 
is not affected or only slightly destroyed during the AGB phase. Thus, the composition of 
PNe is an important tool for understanding the chemical evolution of the Galaxy. Deuterium 
is produced in appreciable quantities only in primordial big bang nucleosynthesis a nd eas- 
ily destroyed in the interiors of stars (jEpstein et al.lll976l ; iTruran fc Cameron! Il97ll ). Many 
studies have been performed through the years in order to determine the current gas-phase 
abundance of deuterium, which can be linked to the primordial abundance of deuterium 



throu gh galactic chemical evolution models (see for e.g. lOliveira et al.ll2006l ; iRomano et al. 



2006 



and references therein). These studies show that, despite infall, the abundance of D 
has decreased over time due to astratio n. The gas-phase D/H is constant throughout the 
Local Bubble (LB), (1.56 ± 0.04) xlO" 5 rtWood et all 120041 ) . Beyond the LB the re is a large 



scatt er in the measurements with values as low a s (0.50 ± 0.16) x 10 5 (6 Car, lAllen et al. 



1992|) and as high as (2.18 ± ai)xl0" 5 {j 2 Vel, ISonneborn et al.ll2000( ). It has been pro- 



posed that the low gas-phase D/H values are due to deuterium being depleted in grains and 
that the tota l present day D/H in the local Galactic disk might be higher than 2.31 xl0~ 5 
( iDraind I2004J : iLinsky et al.ll2006l ). In this work we present a study of D/H along the sight- 
line to WD 0439+466, the central star of the PN Sh 2-216. We argue that the low gas-phase 
D/H ratio derived here, (0.76 ± [j;^) x 10~ 5 , is not due to deuterium being depleted in dust 
grains but is a consequence of the sightline being dominated by astrated, deuterium-poor 
gas associated with the PN. All uncertainties quoted are at the la level. 



2. WD 0439+466 and its Planetary Nebula 



Sh 2-216 is a large and old (3-6xl0 5 vr. iTweedv Napiwotzki 



1992 



Napiwotzki! 1 1999 ) 



low-surface brightness PN at a distance of 129 ±| pc ( [Harris et al.l 120071 . trigonometric 
parallax) in the direction I = 158.49°, b = +0.47° . At this d is tance its apparent size of 
1?6 translates into a linear diameter of ~3.6 pc. iFesen et al.l (Il98ll ) first proposed that 
Sh 2 -216 was a PN b ased on the strength of the [N II] and [S II] emission lines. A few years 
later iReynoldd (119851 ) identified WD 0439+466 (DAO-type white dwarf) as a possible central 
candidate for the Sh 2-216 PN. This identification was difficult because the central star is 
no longer located in the geomet rical cente r of th e nebula, probably due to the interaction of 
the PN with the ambient ISM. IReynoldd (119851 ) derived an expansion velocity < 4 km s _1 
and a temperature T = 9400 ± 1100 K for the nebula. From its distance and proper motion 
Kerber et al.l (120041 ) determined that WD 0439+466 left the geometrical center of the PN 
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Sh 2-216 about 45,000 years ago. 



FUSE and STIS OBSERVATIONS 



WD 0439+466 has been observed by FUSE numerous times through all the apertures 
(programs M105, M107, and P104). We coadded all the calibrated data (Ca lFUSE version 
3.C0) for each aperture individually, following the procedures described in Oliveira et al. 



( 120061 ). The three resulting datasets have exposure times ranging from 10 ksec to 51.5 ksec. 
The FUSE data covers the spectral range 905 - 1185 A with a spectral resolution of ~20 km 



(Moos et al. 


2000; 


Sahnow et al. 


2000) 



The HST/STIS data were obtained with the E140M echelle grating through the 0.2" x 
0.2" aperture and have a resolution of ~5 km s _1 (observations o64d03010 and o64d03020 
with exposure times of 2.4 ksec and 3.0 ksec, respectively). The data, processed with CAL- 
STIS version 2.18, were downloaded from the MAST archive at STScI and no further pro- 
cessing was applied. 



RESULTS AND DISCUSSION 



Column densities of several atomic and molecular species , d etermined from FUSE and 
HST data (following the same procedures of lOliveira et al.ll2006l ). are presented in Table [TJ 



iV(H I) was determined from Lya in conjunction with non-LTE stella r models (with T e & 
= 83200 ± 3300 K, log g = 6. 74 ± 0.19, and log (H e/H) = -1.95 from lNapiwotzki]ll999} ) 
in a manner similar to that of lOliveira et al.l (120061 ) . Our fits yield logiV(H I) = 19.91 ± 
0.06. Figured] (top panel) presents the fit to the Lya region (in blue). The la uncertainties 
of 0.06 dex (in red) combine the statistical uncertainties with the one a ssociated with using 
different stellar models. Our iV(H I) is consistent with that obtained by lRauch et al.l (120061 ) 
using more recent models that also incorporate metals and that lead to T e g = 95,000 K and 
log g = 6.90. (iV(HI) is not very sensitive to the T eS of the stellar model because of the 
large breadth of the Lya profile and the small stellar contribution to the profile.) 

No stellar model was used to determine N(D I). One-absorption-component fits to the 
AA919.1, 920.9 (blended with H 2 J = 6), 922.9 (blended with H 2 J = 3 and HD J = 0), 
and 925.9 (minor blend with H 2 J = 7) D I lines lead to N(D I) = 15.09 ± 0.02. However 



1 The CalFUSE pipeline reference guide is available at 
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the A919.1 D I line is not well fit, possibly due to a blend with an unknown stellar line. 
Removing this line from the fit yields N(D I) = 15.00 ± 0.02. To take into account the 
different fits we adopt for this sightline N(D I) = 15.05 ± 0.07. Figured] (bottom panels) 
presents fits to some of the D I lines used in the analysis. 

At the FUSE and STIS resolutions all the ISM lines display only a single absorption 
component with a common velocity, fhciio ~ 6 km/s; hence in Table [1] we quote total sightline 
column densities. However, the detection of several ionization stages for some of the species 
(e.g. C I, C II, and C IV) indicate that there must be at least three absorption components 
along this line of sight: a cool component traced by molecular hydrogen and C I (amongst 
other species), a photoionized component traced by the high ionization species (C IV, N V, 
and Si IV), and a warm component where the bulk of the species such as C II, O I, N I,S II, 
Fe II, etc. reside. The highly ionized species (no values of iV quoted in Table [T] because 
of blending with the stellar component) likely reside in the photoionized gas associated 
with the planetary nebula where optical emission lines of Ha and O III have b een detected 
rtFesen et al.lll98lh . The LB contributes only log 2V(H) ~ 19.3 Jsfeir et all 1 1993 ) to the total 
neutral hydrogen column density, log iV(H) = 20.30 (or ~10%), so most of the gas sampled 
here is beyond the LB. However, D/H derived below is only an upper limit to D/H in the 
PN because of D and H absorption by foreground gas. 

We derive, for the neutral gas along this sightline, D/H = (0.76 ± Q;Jf)xlO~ 5 , N/H = 
(3.24 ± 0;|)xl0- 5 , and O/H = (0.89 ± [j;^)xl0~ 4 (see Table EJ). Note that this sightline 
contains a significant amount of molecular gas. Therefore, the neutral gas ratios take iV(H 2 ) 
and JV(HD) into account: D/H = (JV(D I) + JV(HD))/(JV(H I) + 2iV(H 2 ) + iV(HD)), etc. 

The D/H ratio is ~2 times smaller than the LB value quoted above but comparable to 
D/H along other sightlines. However, the N/H and O/H rati os, ~2.3 and ~3.9 times smaller 



respectively than the typical IS M ratios of (7.5 ± 0.4) xlO" 5 flMever et al.lll997D and (3.43 ± 
0.15) x 10 -4 (IMeyer et al.lll998l ). are a strong indication that the gas sampled by this sightline 
is different from typical diffuse sightlines. It is unlikely that the O/H ratio is due to depletion 
into dust, as studies of the O/H ratio have found that oxygen is only mildly depleted in hig h 
mean density sightlines, with an average ratio of (2.84 ± 0.12) xl0~ 4 (ICartledge et alj |2004), 
significantly higher than the value derived here. Using equations 14 -18 of PeimbertT(1968 ) 



in conjunction with the optical emission line data from Table 2 of iFesen et al.l (119811 ) we 
estimate, for the ionized gas, O/H ^3.9xl0~ 4 and N/H ~7.3 xl0~ 5 for position C (closest 



to the central star) (using T e = 9400 K from iReynoldsl Il985l ) . These ratios, based on Ha, 
N II, O II, and O III data, are more similar to O /H and N/H ratios typic al of circular PNe 
(~3.6xl0~ 4 and ~1.7xl0~ 4 , respectively, from IStanghellini et al.l 120061 . although a large 
scatter exists). Ionization corrections, not applied here, are then responsible for our low 
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neutral-gas N/H and 0/H ratios but can not explain our D/H ratio, as D I and H I have the 
same ionization fraction due to rapid charge exchange (IBlack fc Dalgarnolll973l ). We note 
that the coupling of O and H i onization fractions i s very strong unless there is appreciable 
ionization of O to higher stages (IJenkins et al.ll2000l ) as is the case along this sightline. Other 
indications that most of the gas along this sightline is associated with the PN come from 
the molecular gas: 1) Molecules such as CO and H2 exist in the envelopes of PNe. These 
molecules are distributed in a clumpy structure and in evolved PNe can survive only inside 
dense clumps or they would be rapidly photoionized and dispersed on time scales shorter 
than the age of the nebula (jTielend Il993l ; iHugginsI l2000l ) . In addition, the fraction of H2 
along this sightline, /h 2 = 0.58 ± 0.06, is not typical of sightlines with such low iV(H) and 



is on e of the highest values in the literature for log iV(H) < 21.0 (see e.g. iGillmon et al. 



20061 . and references therein). Hence, the gas containing the CO and probably most of the 
H 2 is likely in a dense clump in the envelope of the PN. 2) Transitions from high J lines 
(J = 7-9) and HD J = 1 are not typical of diffuse cloud sightlines. These lines are formed 
i n environments with radiation fields much more intense than the average interstellar field 
(ISnow fc McCaHll2006l ). lending support to the idea that the molecular gas is close to the 
central star of the PN. 

Although depletion of D into dust grains cann ot be ruled out we think it is unlikely be- 
cause: 1) Sh 2-216 is much younger (~ 0.3-0.6 Myr. iTweedy &: Napiwotzkilll992l ; lNapiwotzki 



19991 ) than th e timescale fo r deuterium depletion into dust grains (~2.3 Myr for D/H to de 
crease by 11 e, lDrainell2006l ) . 2) Extreme enrichment of the grains relative to the gas, requires 
^dust < 90 K in order to explain the lowest D/H ratios observed (Drainel 2006 ; Linsky et al. 



20061 ) but dust in PNe is typically warmer (90-150 K. Ivan Hoof et alJbood ) 



The low D/H, N/H, and O/H ratios combined with the detection of CO, high H 2 and 
HD J lines, and the high /h 2 are consistent with the idea that the gas along this sightline is 
dominated by gas that is associated with the PN. In this scenario the low D/H ratio is not 
a consequence of depletion into dust grains (although this cannot be completely ruled out) 
but of the gas in the nebula being deuterium-poor as a result of astration. 

Even though astration has reduced the primordial abundance of deuterium to the lower 
values found today in the Galaxy, this would be the first time that we have seen directly the 
effects of astration in the D/H ratio of any single sightline (see e.g. iLinsky et al.ll2006l . and 
references therein, for a summary of D/H ratios and pssobile explanations for the low values). 
This work raises the possibility that other low D/H ratios could be explained sometimes by 
the presence of astrated gas that has not yet mixed with the ambient gas, particularly for 
nearby sightlines where this effect is less likely to be averaged out. 
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Table 1. Adopted Column Densities (cm 2 ) 





Lop- N 




Species 


UUg, 1 V 


H I 


19.91 ± 0.06 


Fe II 




14.11 ± 0.06 


D I 


15.05 ± 0.07 


H 2 


(J = 


0) 


19.54 ± 0.05 


C I 


14.60 ± 0.10 


H 2 


(J = 


1) 


19.37 ± 0.05 


CI* 


14.40 ± 0.06 


H 2 


(J = 


2) 


18.6: 


CI** 


14.02 ± 0.08 


H 2 


(J = 


3) 


18.0: 


C II 


>14.40 


H 2 


(J = 


4) 


16.7: 


CII* 


>14.25 


H 2 


(J = 


5) 


16.5: 


NI 


15.81 ± 0.07 


H 2 


(J = 


6) 


15.00 ± 0.25 


Nil 


>14.17 


H 2 


(J = 


7) 


14.85 ± 0.08 


01 


16.25 ± 0.05 


H 2 


(J = 


8) 


13.70 ± 0.10 


Al II 


>12.58 


H 2 


(J = 


9) 


13.30 ± 0.05 


Si II 


>14.08 


H 2 






19.76 ± 0.04 


Si III 


>12.67 


H 






20.30 ± 0.03 


p II 


>12.70 


HD (J = 


= 0) 


14.54 ± 0.10 


S I 


13.45 ± 0.08 


HD (J = 


= 1) 


13.72 ± 0.05 


S II 


>14.90 


HD 






14.60 ± 0.09 


CII 


>12.71 


D = 


-- D I + HD 


15.18 ± 0.06 


Ar I 


>13.53 


CO 






13.78 ± 0.07 



Note. — Quoted uncertainties are la. For H 2 J = 2-5 
the column densities are uncertain; these are flagged with 
':'. We use 7V(H 2 ) ~ N(J = 0) + N ( J = 1) since these 
levels contain most of the H 2 along the sightline. iV(H) = 
iV(H I) + 2iV(H 2 ) + iV(HD) 
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Table 2. Ratios of Column Densities 



Species Ratio Species Ratio 



D/H 


(0.76 ± 


8:i?)xio- 


-5 


N/O 


0.36 ± °;« 


DI/HI 


(1.38 ± 


8:2?) x io- 


-5 


HD/2H 2 


(3.39 ± 8'6?)xl0- 6 


O/H 


(0.89 ± 


°o.\l) x 10" 


-4 


HD/D 


0.26 ± ™ 7 


N/H 


(3.24 ± 


8:1) x io- 


-5 


CO/H 2 


(i.05±g:S)xio- 8 


Fe/H 


(0.65 ± 


8:£) x 10" 


-6 


/h 2 


0.58 ± 0.06 


D/O 


(8.51 ± 


1:11) x io- 


-2 


(n H ) (cm" 3 ) 


0.54 



Note. — Quoted uncertainties are la. N(D) = N(D I) + iV(HD), 
N(R) = N(R I) + 2xA r (H 2 ) + AT(HD), and / Ha = 2x AT(H 2 )/7V(H). 
For the ratios above we assume jV(O) « iV(0 I), iV(N) « iV(N I), 
and iV(Fe) « iV(Fe II). 
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Fig. 1. — Top Panel: Fit to H I Lya (bl ue), obtained with a model with T e fj = 83,200 K, 
log g = 6.74, and log (He/H) = -1.95 JNapiwotzkil Il999h . yielding log iV(H I) = 19.91. 
1 cr uncertainties of 0.06 dex are represented by red. Bottom Panels: Fits to some of the 
D I lines used in the analysis (most and least blended D I lines used, left and right panels, 
respectively). H I is represented by a black dashed line, D I in red, O I in green, Fe II in 
magenta, H 2 in blue, HD in orange, and CO in dashed grey. The overall fit is represented 
by a thick grey line. The FUSE channel used in each plot is given at the bottom left. S2AL 
corresponds to LWRS SiC 2A. 



